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[i]  The  bubble  cloud  depth  and  its  correlation  with  extreme 
winds  are  key  elements  of  bubble-mediated  gas  injection, 
which  are  critical  to  the  determination  of  the  global  gas 
budgets.  The  characteristics  of  bubble  cloud  depth  were 
examined  from  measurements  collected  during  the  passage  of 
a  category-4  hurricane  with  winds  up  to  50  m  s  1 .  The  bubble 
cloud  depth  increases  linearly  with  wind  speed  for  winds  less 
than  35  m  s  ’.  Our  findings  are  consistent  with  previous 
observations  at  low  to  moderate  wind  speeds.  However, 
the  rate  of  increase  is  reduced  significantly  at  winds  higher 
than  35  m  s  1 .  Citation:  Wang,  D.  W.,  H.  W.  Wijesekera,  W.  J. 
Teague,  W.  E.  Rogers,  and  E.  Jarosz  (2011),  Bubble  cloud  depth 
under  a  hurricane,  Geophys.  Res.  Lett.,  38,  LI 4604,  doi:  10. 1029/ 
201 1GL047966. 

1.  Introduction 

[2]  Bubble  clouds  generated  by  breaking  waves  play  an 
important  role  in  gas  exchange  between  the  atmosphere  and 
upper-ocean.  Formation  of  bubble  clouds  during  storms  is 
responsible  for  injecting  nitrogen  and  oxygen  gases  that  can 
supersaturate  the  ocean  [Vagle  et  al.,  2010].  Field  data  from 
Hurricane  Frances  (2004)  indicated  that  the  complete  bubble 
injection  process  at  high  winds  becomes  an  important  factor 
in  the  air-sea  gas  exchange  process  [D  'Asaro  and  McNeil, 
2007;  McNeil  and  D  ’Asaro,  2007].  Recently,  oxygen  and 
nitrogen  gas  injection  coefficients  are  shown  to  be  related  to 
the  bubble  cloud  depth  during  winter  storms  over  the  north¬ 
west  Pacific  Ocean  [Vagle  et  al.,  2010],  Observations  of 
bubble  cloud  depth  evolution  directly  under  hurricane  paths 
are  extremely  rare  due  to  the  logistics  of  making  such  mea¬ 
surements.  As  a  result,  characteristics  of  bubble  cloud  depth 
and  their  correlations  with  wind  forces  are  not  well  under¬ 
stood.  These  are  key  elements  to  the  understanding  of  the 
injection  of  weakly  soluble  gases  such  as  nitrogen  and  oxy¬ 
gen,  which  is  becoming  more  critical  to  the  determination  of 
the  global  gas  budgets,  and  ultimately  to  the  global  climate 
since  the  number  and  intensity  of  storms  appears  to  be 
increasing  over  the  past  and  future  decades  [Emanuel,  2005; 
Bender  et  al.,  2010]. 

[3]  A  unique  data  set  of  acoustic  backscatter  was  collected 
under  a  category-4  hurricane,  Ivan  in  2004,  on  the  continental 
shelf  in  the  northeastern  Gulf  of  Mexico  as  part  of  the  Naval 
Research  Laboratory’s  (NRL’s)  Slope  to  Shelf  Energetics  and 
Exchange  Dynamics  (SEED)  project  [Teague  et  al.,  2007]. 
Hurricane  Ivan  passed  directly  over  the  SEED  mooring  array 
on  around  0000  UTC,  16  September  2004  (Figure  1)  [Teague 
et  al.,  2007].  Although  the  main  observations  from  SEED  are 
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current  velocity  profiles,  the  data  set  also  contains  profiles 
of  acoustic  backscatter  from  six  bottom-mounted  acoustic 
Doppler  current  profilers  (ADCPs)  near  29.4°N,  88°W.  This 
study  examines  acoustic  backscattering  intensity  as  a  mea¬ 
sure  of  bubble  concentration  along  the  path  of  Hurricane  Ivan. 
The  main  objective  of  this  study  is  to  quantify  the  bubble 
cloud  depths  and  their  relationships  to  wind  speed. 

2.  Measurements  and  Analysis 

[4]  During  NRL’s  SEED  project,  six  Bamy  moorings  were 
deployed  on  the  outer  continental  shelf  along  60-m  and  90-m 
isobaths.  Ml,  M2,  and  M3  were  on  the  60-m  isobath  and  M4, 
M5,  and  M6  were  on  the  90-m  isobath  (Figure  1).  Moorings 
were  spaced  approximately  15  km  apart  on  each  isobaths. 
Each  mooring  carried  a  300-kHz  RD  Instruments  Workhorse 
ADCP  and  a  wave/tide  gauge.  ADCPs  rested  about  0.5  m 
above  the  ocean  floor  with  four  transducers,  each  with  an 
incidence  angle  of  20  degrees.  Echo  intensity  profiles  with 
2-m  cell  resolution  were  ensemble  averages  of  55  pings  that 
were  transmitted  over  a  15-min  period  for  each  transducer 
beam.  The  depth  reference  for  ADCP  profiles  was  taken  as 
the  distance  from  the  mean  sea  surface,  which  was  estimated 
from  the  15-minute  averaged  bottom  pressure  with  a  cor¬ 
rection  of  sea  surface  pressure  [Wang  et  al.,  2005].  The  echo 
intensity  was  converted  to  volumetric  backscatter  strength 
Mv  in  m~'  by  the  use  of  the  sonar  equation  with  calibration 
coefficients  provided  by  the  manufacturer,  Teledyne  RD 
Instruments  (http://www.rdinstruments.com/).  Details  of  the 
backscatter  equation  and  coefficients  are  given  by  Deities 
[1999].  The  backscattering  strength  is  also  expressed  as 
Sy  =  1 0  log i o(/V/v),  where  Sv  is  in  dB.  The  backscattering  data 
from  the  near-surface  ADCP-cells  were  likely  contaminated 
by  the  side-lobe  effect,  and  therefore  we  excluded  data  from 
the  3.6-m  and  5.4-m  cell  depths  for  the  60-m  and  90-m 
deployments,  respectively.  The  cells  near  the  mean  sea  level 
would  be  above  or  below  of  the  ocean  wave  surface  during 
data  acquisition.  In  this  study,  we  excluded  backscattering 
data  from  the  cells  that  were  within  the  amplitude  of  the 
surface  waves  relative  to  the  mean  sea  level.  Here,  the  wave 
amplitude  was  taken  as  the  half  of  significant  wave  height 
(Hsl 2).  In  addition,  the  averaged  ADCP  backscatter  profiles 
were  based  on  a  depth  reference  with  respect  to  mean  sea 
level,  which  led  to  a  larger  backscatter  Sv  than  those  based  on 
a  depth  reference  with  respect  to  the  wave  surface  [Thorpe, 
1982], 

[5]  Backscatter  profiles  for  each  SEED  mooring  were 
constructed  by  averaging  of  the  profiles  of  four  transducer 
beams.  The  averaged  profile  of  acoustic  backscatter  Mv 
shows  an  exponential  decrease  with  depth,  and  can  be 
expressed  as  My  =  M0  exp(-z/Z), where  z  is  the  depth  below 
mean  sea  level;  M0  is  the  backscatter  level  at  z  =  0;  and  L  is  the 
characteristic  (“e-folding”)  length  scale,  which  detenuines 
the  rate  of  decrease  of  Mv  with  increasing  depth  [Thorpe, 
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Figure  1.  Path  of  Hurricane  Ivan  (red  dashed  line).  Mooring  locations  at  60  m  and  90  m  isobaths  are  marked  in  solid  circles 
(M1-M6).  The  yellow  triangle  is  the  NDBC  buoy  42040.  The  insert  shows  a  BARNY  mooring  before  deployment  (lower  left). 
Contour  depths  are  in  meters. 


1982,  1984;  Trevorrow,  2003;  Vagle  et  al.,  2010].  The 
entrainment  depth  of  bubbles  can  be  estimated  from  a  given 
backscatter  profile  by  knowing  the  noise  level  of  the  mea¬ 
surement.  Here  the  bubble  cloud  depth  Zh  was  determined  to 
represent  the  lower  boundary  of  the  bubble  layer,  where  the 
backscatter  intensity  exceeds  an  empirically-determined 
noise  threshold  MT.  In  this  study,  MT  =  10  6  (or  Sv  =  -60  dB) 
which  was  the  noise  threshold  prior  to  the  arrival  of  Ivan. 
Similar  techniques  have  been  used  to  determine  bubble 
cloud  depths  with  slightly  different  threshold  levels,  -50dB 
[Trevorrow,  2003]  and  -70dB  [Vagle  et  al,  2010],  The 
e-folding  length  scale,  L  was  estimated  by  fitting  the  expo¬ 
nential  profile  between  Zh  and  the  acceptable  shallowest 
depth. 

[6]  Ivan  approached  the  mooring  site  as  a  category-4  hur¬ 
ricane  with  a  steady  speed,  direction,  and  intensity.  Its  eye 
passed  through  the  array  and  hence  its  extreme  winds  were 
prevalent  over  all  six  moorings  [Teague  et  al.,  2007].  The 
diameter  of  the  eye  was  about  50  km.  The  center  of  the  eye 
passed  directly  over  the  M2  and  M5  moorings,  while  M3 
and  M6  were  located  slightly  to  the  right  side  of  the  hurricane, 
and  Ml  and  M4  were  located  slightly  to  the  left  side  of  the 
hurricane  (Figure  1).  Ivan  produced  extreme  waves  with 
significant  wave  heights  higher  than  1 5  m  [  Wang  et  al. ,  2005] 
and  currents  of  about  2ms1  [Teague  et  al. ,  2007 ;  Wijesekera 
et  al.  2010].  During  the  passage  of  Ivan  on  September  16, 
2004,  the  water  above  each  SEED  mooring  (Figure  1)  was 
well  mixed  according  to  the  small  water  temperature  differ¬ 
ence  (<2°C)  between  the  bottom  temperature  from  Bamys 
and  the  sea  surface  temperature  from  satellite  and  National 
Data  Buoy  Center  (NDBC)  buoy  42040  [e.g.,  Teague  et  al., 
2007].  The  surface  wind  over  the  mooring  site  was  recon¬ 
structed  by  combining  wind  data  from  the  NDBC  buoy  42040 
and  model  winds  from  the  post-storm  wind  analysis  by  the 
Hurricane  Research  Centcr(HRD)  [Powell  et  al.,  1998;  Wang 
et  al.,  2005;  Ardhuin  et  al.,  2010].  Wave  fields  at  different 


mooring  locations  were  obtained  from  the  wave  model 
“Wave Watch  III”  [Ardhuin  etal,  201 0]  based  on  HRD  winds 
and  NDBC  wave  buoy  data. 

[7]  As  seen  in  Figure  2a,  the  hurricane  passage  over  M5  can 
be  divided  into  four  stages.  They  are:  (I)  the  arrival  of  the  eye- 
wall  with  increasing  wind,  (II)  the  arrival  of  the  eye  with 
decreasing  wind,  (III)  the  departure  of  the  eye  with  increasing 
wind  until  the  arrival  of  backside  of  eye-wall,  and  (IV)  the 
departure  of  the  eye-wall  with  decreasing  wind.  Hurricane 
Ivan’s  eye  took  about  four  hours  to  cross  the  mooring  array. 
Wind  speeds  decelerated  from  45  m  s  1  to  15  m  s  1  and  then 
accelerated  back  to  about  40  m  s  1 .  The  backscatter  obser¬ 
vations  from  all  six  moorings  were  from  Ivan’s  front-right 
and  rear-left  quadrants.  Figure  2b  shows  relative  locations  of 
M5  to  the  eye  of  hurricane  and  the  forward  direction. 

[8]  A  dramatic  change  in  the  ADCP  backscattering 
strength  was  found  during  the  passage  of  Ivan.  As  shown  in 
Figure  2a,  prior  to  the  arrival  of  Ivan,  intense  backscatter 
induced  by  bubble  clouds  was  limited  to  about  10  m  beneath 
the  mean  sea  surface.  The  bubble  cloud  depth  increased  with 
increasing  wind  speed  and  wave  height.  As  a  result,  a  back¬ 
scatter  intensity  of  about  -70dB  can  be  found  as  deep  as  40  m 
when  Ivan’s  eye-wall  reached  M5  with  peak  winds  of  about 
45  m  s'1.  The  bubble-induced  backscatter  layer  decreased 
and  increased  rapidly  when  the  eye  crossed  M5  (stages  II 
and  III),  during  which  the  ocean  was  experiencing  a  rapid 
decrease  and  increase  of  winds,  surface  waves,  and  near¬ 
surface  currents.  The  bubble  cloud  layer  increased  rapidly  as 
the  eye-wall  of  Ivan  crossed  the  mooring  site  for  the  second 
time.  Similar  relationships  between  the  layer  of  high  back¬ 
scatter  intensity  and  local  winds  were  observed  at  the  rest  of 
the  mooring  sites.  Three  selected  ADCP  backscatter  profiles 
from  M5  at  wind  speeds  of  19.9  m  s_1,  29.6  m  s_1  and  41  m 
s_1  are  shown  in  Figure  2c.  The  bubble  cloud  depth  Zb  is 
determined  as  the  depth  at  which  Sv  =  -60dB.  The  backscatter 
strength  between  the  shallowest  uncontaminated  depth  and 
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Figure  2.  (a)  (top)  Surface  wind  at  1 0  m  height  (multi-  colored  solid  line)  and  significant  wave  height  Hs  (blue  dashed  line)  at 
M5  during  Ivan’s  passage.  Red,  black,  green,  and  magenta  wind  segments  represent  the  four  stages  of  Ivan’s  passage,  respec¬ 
tively.  (bottom)  Contour  plot  of  ADCP  backscatter  profile  from  M5  during  Ivan’s  passage.  -60dB  contour  line  is  highlighted 
as  magenta  line,  (b)  Polar  plot  shows  the  relative  location  of  M5  with  respect  to  the  Ivan’s  center  and  its  forward  propagation 
direction  (0  degree)  at  the  times  of  backscatter  data  acquired  within  200  km  of  Ivan  center  (after  1630  UTC  15  September). 
Radial  distance  is  in  km.  Color  symbols  represent  the  four  stages  as  defined  in  Figure  2a.  (c)  Three  ADCP  backscatter  depth 
profiles  from  M5  at  12:45  UTC  (green),  19:45  UTC  (cyan),  and  22:30  UTC  (magenta),  September  15  are  selected.  Winds  are 
19.9ms  1,29.6ms  1  and41  ms  ’,  and  significant  wave  heights  are  8.1  m,  14.6mand  13.8  m,  respectively.  The  near-surface 
sections  which  were  disregarded  in  the  computation  are  not  shown.  Black  dotted  line  represents  the  fitted  exponential  decrease 
with  depth  Mv  =  M0  exp(-z JL)  for  each  profile,  (d)  Bubble  cloud  depth  Zb  vs.  wind  speed  for  the  four  stages  of  Ivan’s  passage 
from  M5.  The  color  symbols  represent  the  four  stages  as  defined  in  Figure  2a. 


Zb  decays  exponentially  with  depth  and  is  well  represented  by 
Mv  =  Mq  exp(-z JL). 

[9]  The  scatter  plot  ofZb  versus  wind  is  shown  in  Figure  2d. 
For  winds  less  than  15  m  s-1,  Zb  is  relatively  independent 
of  wind  due  to  the  combined  effects  of  a  coarse  vertical 
cell  size  (2  m)  and  the  side-lobe  effect  (5.4  m)  of  the 
ADCP.  Therefore,  we  focus  our  analyses  for  winds  larger 
than  15  ms  1 .  During  the  arrival  of  the  eye-wall  (stage  I), 
Zb  increased  with  increasing  wind  speed.  During  stages  II 
and  III,  the  bubble  cloud  depth  first  decreased  and  then 
increased  as  a  result  of  the  dramatic  variation  of  wind  speeds 
across  the  eye.  The  bubble  cloud  depth  is  closely  related  to 
the  winds,  but  the  rate  of  change  in  bubble  cloud  depth  was 
generally  faster  inside  the  eye-wall  than  the  eye-wall  arrival 
stage  (stage  I).  The  differences  in  the  correlation  of  Zb  with 
winds  at  the  four  hurricane  passage  stages  could  be  due 
to  variations  of  breaking-wave  characteristics.  The  aver¬ 
aged  relationship  between  Zb  and  wind  speed  during  the 
growing  stage  (i.e.,  stage  I,  Figure  2d)  of  the  hurricane  was 
examined  by  averaging  Zh  into  wind  speed  bins  using  data 
from  all  six  moorings.  The  bin-averaged  bubble  depth  (Zb) 
versus  wind  speed  is  shown  in  Figure  3a.  (Zb)  is  closely  related 
to  wind  speed  and  increased  from  about  7.5  m  to  about  2 1  m  as 
wind  speed  increased  from  15  m  s  1  to  35  m  s  The  observed 
relationship  is  consistent  with  results  from  Vagle  et  al.  [2010] 
based  on  backscatter  profiles  observed  for  winds  less  than 
25  m  s  1  (dashed  blue  line  in  Figure  3a).  For  winds  higher 
than  35  m  s  \  the  rate  of  change  of  (Zb)  became  smaller  in 
spite  of  increasing  winds.  This  coincides  with  the  reduction  of 
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Figure  3.  Bin-averaged  (a)  bubble  cloud  depth  (Zb)  and 
(b)  e-folding  length  ( L )  versus  wind  speed.  Error  bar  lines 
indicate  one  standard  deviation.  The  thick  dash  lines  in 
Figures  3a  and  3b  represent  the  linear  relationship  from 
Vagle  et  al.  [2010], 
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the  surface  drag  coefficient  at  high  winds  as  reported  by 
Powell  et  al.  [2003]  and  Jarosz  et  al.  [2007]. 

[10]  The  e-folding  length  L  computed  from  the  backscatter 
profile  reflects  the  rate  of  the  decrease  in  bubble  concentra¬ 
tion  with  an  increasing  depth.  The  e-folding  length  appears  to 
be  related  to  the  characteristic  length  scale  associated  with  the 
turbulent  diffusion  coefficient  in  the  surface  layer  [ Thorpe , 
1984].  As  shown  in  Figure  3b,  the  bin-averaged  e-folding 
length  (L)  increased  from  about  0.8  m  at  winds  of  1 5  m  s  1  to 
about  2.8  m  at  winds  of  35  m  s  1 .  The  observed  growth  rate  of 
(L)  is  qualitatively  consistent  with  the  linear  relation  from 
Vagle  et  al.  [2010]  based  on  their  measurements  for  winds 
under  25ms1.  Similar  to  (Zb) ,  (L)  flattened  for  winds  higher 
than  35  m  s  \  Estimated  e-folding  lengths  from  Vagle  et  al. 
[2010]  are  higher  than  our  estimates  by  a  factor  of  1.5-2. 
These  differences  could  be  because  the  backscatter  profiles 
from  Vagle  et  al.  [20 1 0]  were  based  on  a  depth  reference  with 
respect  to  the  wave  surface,  which  can  result  in  a  higher 
e-folding  length  scale  than  a  depth  reference  based  on  a  mean 
sea  level  as  used  here  [e.g.,  Thorpe,  1982]. 

3.  Discussion  and  Summary 

[n]  We  examined  and  quantified  the  evolution  of  bubble 
cloud  depth  based  on  measurements  during  the  passage  of  a 
category-4  hurricane  with  winds  up  to  50  m  s  .  We  found 
that  the  bubble  cloud  depth  and  e-folding  length  scale 
increase  with  wind  speed.  However  the  growth  rates  of  (Zb) 
and  ( L )  were  much  reduced  when  winds  exceeded  35  m  s_1. 
The  flattening  of  (Zb)  at  winds  higher  than  35  m  s  1  may 
explain  why  gas  injection  flux  peaked  at  winds  of  about  37  m 
s  1  during  observations  of  gas  transfer  rates  by  D  ’Asaro  and 
McNeil  [2007]  from  Hurricane  Frances.  These  findings  on 
the  wind  speed-dependent  (Zb)  and  ( L )  lead  us  to  postulate: 
(a)  wind  driven  turbulences  increase  Zb,  (b)  at  high  winds,  Zb 
and  L  are  related  to  surface  drag  and  (c)  the  reduced  rate  of 
bubble  cloud  depth  and  e-folding  length  at  high  winds  could 
be  due  to  the  combination  of  the  reduction  of  surface  drag 
[Jarosz  et  al,  2007;  Powell  et  al.,  2003]  and  the  complete 
bubble  dissolving  under  hydrostatic  pressure  [D  Asaro  and 
McNeil,  2007;  McNeil  and  D  Asaro,  2007].  This  also  implies 
that  the  parameterization  of  air-sea  gas  exchange  based  on 
wind  speed  alone  should  be  limited  to  winds  up  to  about 
3  5  m  s~ 1 .  Extending  the  wind-speed  dependent  relation  of  gas 
exchange  to  high  winds  could  result  into  an  overestimation. 
The  dynamics  of  air-sea  momentum  transport  (such  as  the 
complexity  of  wind-stress  dependence  on  the  surface  drag 
coefficient)  at  the  sea  surface  must  be  considered  for 
parameterization  of  gas  exchange  at  high  winds. 
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